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Abstract 

Fragmentation functions for charged particles in — > qq events have been measured 
for bottom (b), charm (c) and Hght (uds) quarks as well as for all flavours together. The 
results are based on data recorded between 1990 and 1995 using the OPAL detector at 
LEP. Event samples with different flavour compositions were formed using reconstructed 
D*^ mesons and secondary vertices. The = ln(l/xp) distributions and the position 
of their maxima Eire also presented separately for uds, c and b quark events. The 
fragmentation function for b quarks is significantly softer than for uds quarks. 
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1 Introduction 



Experimental measurements of the inclusive momentum distribution of charged particles in 
e"'"e~ collisions provide important insight into the process of how quarks turn into hadrons. 
This distribution is commonly normalised to the total hadronic cross section cr*°* and presented 
as a function of the scaled momenta Xp = 2ph/\/s of the charged hadrons, where ^/s is the 
centre-of-mass energy. In this form, the spectrum is usually referred to as the fragmentation 
function and can be obtained experimentally from the total number of hadronic final states, 
iVevent, and the number of charged particles in each Xp bin, A^track(a^p): 



The charged particle momentum spectrum can also be studied as the distribution of ^p = 
ln(l/a;p). The C,p distribution emphasises the low momentum component and the Xp distribu- 
tion the high momentum component of the momentum spectrum. 

In the naive quark parton model, the scaled momentum distribution is expected to be 
independent from the centre-of-mass energy. A violation of this scaling is expected due to 
gluon radiation in the final state. Experimentally, scaling violation in fragmentation functions 
had indeed been observed by combining measurements at different centre-of-mass energies and 
could be used to determine as [jl]]. The position of the maximum of the ^p distribution, ^o, has 
been studied in the past in various experiments (see for example 0] and references therein). 
The energy dependence of the position of the maximum provides an important test of the 
QCD prediction for the emission of soft gluons 0. 

In events with a heavy primary quark, the possibility of cascade decays of bottom or charm 
hadrons results in more particles sharing the same energy than in light quark events and a 
softer momentum spectrum can be expected. Since the flavour composition of the primary 
quarks in e"'"e~ — ^ qq is predicted by the electroweak theory to change with centre-of-mass 
energy, this flavour dependence of the momentum spectra affects the energy dependence of the 
Xp and distributions of the inclusive event sample. To correct for this contribution, in |Q 
not only the inclusive fragmentation function was measured but also fragmentation functions 
in event samples with different flavour compositions were studied. In measurements of 
fragmentation functions in samples with different flavour composition were used to extract 
flavour dependent fragmentation functions for events with primary light (uds), charm (c) or 
bottom (b) quarks. 

Here we present a measurement of flavour dependent fragmentation functions, based on the 
methods developed for the OPAL measurement of charged particle multiplicities in uds, c and 
b quark events ^. Events were divided into two hemispheres by a plane perpendicular to the 
thrust axis. Secondary vertices and reconstructed D*^ mesons were used to tag hemispheres 
to create samples of events with different quark flavour mixtures (Section 2). To reduce the 
biases induced by the tagging, the measurement of the fragmentation functions was based on 
the momentum spectrum of charged particles in the event hemisphere opposite to the tag. 
Corrections for hemisphere correlations and for distortions due to detector effects are described 
in Section 3. The flavour dependent Xp and C,p distributions were obtained from a simultaneous 




1 da^ 



1 A^track('''p) 



(1) 



dxp 
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fit to the momentum spectra of tlie different hemisphere samples (Section 4). For the first 
time, the measured position of the maximum of the distribution, ^o, is presented separately 
for uds, c and b events. A measurement of the inclusive distribution of all five flavours was 
also performed, based on the track momentum spectrum of all events, i.e., without considering 
any flavour tagging. 



2 Selection and event tagging 

A complete description of the OPAL detector can be found elsewhere p[. This analysis 
relied on the precise reconstruction of charged particles in the central detector, consisting of a 
silicon microvertex detector, a vertex drift chamber, a large volume jet chamber and chambers 
measuring the z-coordinateQ of tracks as they leave the jet chamber. 

This analysis used data recorded with the OPAL detector in the years 1990 to 1995 
at centre-of-mass energies around 91.2 GeV comprising an integrated luminosity of about 
177 pb^^. Z° decays were selected using the criteria described in 0. To ensure that most 
charged particles were well contained in the detector, the polar angle of the thrust axis was 
required to satisfy | cos 6'thrust I < 0.8. To reduce systematic errors in the application of the 
secondary vertex tag, it was only based on a homogeneous data sample taken in the year 1994, 
representing an integrated luminosity of about 34 pb~^. The full integrated luminosity was 
used in the case of the D*^ meson tag. 

Charged tracks used in the measurement of the fragmentation function were required to 
have a measured momentum in the x-y plane, pt, of at least 0.150 GeV/c and to satisfy 
\do\ < 0.5 cm, where do is the distance of closest approach to the origin in the x-y plane. 

Simulated hadronic decays were generated with the Jetset 7.4 Monte Carlo program 
T0| tuned to OPAL data [^. The events were passed through a detailed simulation of the 



OPAL detector [|T2| and processed using the same reconstruction and selection algorithms as 
the data. 



2.1 Secondary vertex tag 

Samples with varying purity of b quark events were selected by reconstructing secondary 
vertices, following the procedure described in 0. Events were divided into two hemispheres 
by the plane perpendicular to the thrust axis and comprising the interaction point. Jets 
were reconstructed by combining charged tracks and electromagnetic clusters not associated 
to tracks, using the scaled invariant mass algorithm described in [|13] with the JADE-EO 
recombination scheme and the invariant mass cut-off being set to 7GeV/c^. A vertex fit was 
then attempted in the highest energy jet in each hemisphere separately. Each track used in 
these vertex fits was required to have at least one hit in the silicon microvertex detector. All 

^The OPAL coordinate system is defined with positive z along the electron beam direction and with positive 
X pointing towards the centre of the LEP ring. The polar angle 6 is defined relative to the +z axis and the 
azimuthal angle (j) relative to the +x axis. 
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such tracks in the jet were fitted to a common vertex point in the x-y plane and the track with 
the largest contribution to the was removed if this contribution was greater than four. The 
remaining tracks were then refitted until either all tracks contributed less than four to the 
or there were fewer than four remaining tracks. For each successfully reconstructed secondary 
vertex, the projected decay length L in the x-y plane with respect to the primary vertex was 
calculated, where the primary vertex was reconstructed from all tracks in the event together 
with a constraint to the average beamspot position as in [|14|. 

The decay length significance, i.e. the decay length divided by its uncertainty L/ctl was 
used to obtain three event samples k of varying b flavour purity. According to the simulation, 
the b purities in these samples vary from 11% to 90% (see Table Fig [l|(a) shows the 
distribution of the decay length signiflcance in data and in Monte Carlo. 





-10.0 < L/aL < 1.0 


1.0 < L/aL < 5.0 


5.0 < L/aL < 50.0 


Number of hemispheres 


940 275 


268 500 


117665 


uds quark fraction fljf^ 


0.71 


0.39 


0.04 


c quark fraction 


0.18 


0.23 


0.06 


b quark fraction 


0.11 


0.38 


0.90 



Table 1: Number of tagged hemispheres in data and the flavour composition derived from the 
Monte Carlo simulation in three decay length significance regions. 



2.2 D*^ meson tag 

Event samples with an enriched c quark contribution were obtained by reconstructing D*"*^ me- 
son candidates. D*^ candidates were selected via the decayQ D*+ K^7r+7r+ closely following 
the procedure described in 0: 

• A subset of tracks was selected that have pt > 0.250 GeV/c and |(io| < 0.5 cm. 

• Candidates of D° — > K~it~^ decays were selected by taking all combinations of two 
oppositely charged tracks, with one of them assumed to be a pion and the other assumed 
to be a kaon. D*^ candidates were selected by combining D° candidates with a third 
track. This 'slow pion' track was required to have the same charge as the track presumed 
to be the pion in the D" decay. 

• The probability that the measured rate of energy loss, dE/dx, for the kaon candidate 
track was consistent with that expected for a real kaon was required to be greater than 
10%. 

• At least two of the three tracks were required to have either z-chamber hits or a polar 
angle measurement derived from the point at which the track has left the jet chamber. 

^Throughout this paper, charge conjugate particles and decay modes are always implied. 
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• The invariant mass of the D° candidate was required to be between 1.790 GeV/c^ and 
1.940 GeV/c^ and the mass difference between the D° and the D*^ candidate, AM, was 
required to be between 0.142 GeV/c^ and 0.149 GeV/c^. 

• Making use of the fact that real D*^ mesons decay isotropically in their rest frames 
whereas combinatorial background is peaked in the forward and the backward direction, 
the following cuts were applied: |cos6'*| < 0.8 for xd* < 0.5 and |cos6'*| < 0.9 for 
Xn* > 0.5, where 6* is the angle between the kaon in the D° rest frame and the direction 
of the D° in the laboratory frame and xd* is the scaled energy of the D*^, i.e., xo* = 
2Ed.±/v^. 



To provide samples with differing charm purity, the data were divided into three xd* regions. 
To evaluate their flavour composition, AM distributions obtained without the cut on AM 
were fitted with a Gaussian for the signal and a function Aexp(— i?AM)(AM/m^ — l)*" for 
the background [|15]. The signals, together with the fitted functions are shown in Fig. |I](b)-(d). 
These fits were used to determine the fraction of background, /^*^ in each D*^ sample in the 
signal AM region. The results of these fits are summarised in Table ^. 





0.2 < Xd* < 0.4 


0.4 < Xd* < 0.6 


0.6 < Xd* < 1.0 


Number of D*^ candidates 
combinatorial background fraction f^'^ 
c quark fraction 
b quark fraction 


5109 

0.57 
0.22 ±0.06 
0.78 ±0.06 


1951 

0.36 
0.50 ±0.06 
0.50 ±0.06 


985 

0.24 
0.90 ±0.04 
0.10 ±0.04 



Table 2: Number of D*^ candidates, the fitted background fraction and tlie flavour composition 
of events witli a genuine D*"^ as taken from jl^/ in three xd* regions. 



The selected samples of D*^ candidates have three components: genuine D*^ mesons from 
b quark decays, genuine D*^ mesons from c quark decays and combinatorial background. No 
other sources of D*^ candidates were considered since Monte Carlo simulations predicts that 
only 0.3% of D*^ mesons with xd* > 0.2 are produced via gluon splitting in light quark events 
ll5| . To evaluate the effect of the contribution from fake D*^, a side-band sample was se- 
lected by requiring that the two pions of the D*^ candidates had opposite charge and that 
0.150 GeV/c^ < AM < 0.170 GeV/c^. Once this contribution was taken into account, the 
flavour composition of the D*^ samples was taken as the fractions and V]^ of genuine 



D* mesons originating from a primary c quark and b quark as measured in |T^. While in 



TSi the fractions and were derived for D*^ candidates after corrections for detector effi- 



ciency and acceptance were made, they were applied in this analysis to uncorrected data. No 
modifications were made since tests with Monte Carlo simulated events showed no significant 
fiavour dependence of these corrections. 
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3 Corrections 



The track momentum distributions in the hemispheres opposite to the secondary vertex or 
D*^ tag were measured. Six distributions (label k) were obtained, corresponding to the three 
decay length regions and the three xd* regions. To obtain fragmentation functions from these 
distributions, three sets of corrections were applied. Firstly, a correction was made to take 
into account track momentum resolution and reconstruction efficiency. Secondly, the effects 
due to the event selection and the correlation between hemispheres were accounted for. In 
addition, the measured track momentum spectra in the D*^ tag samples were corrected for 
the contribution of fake D*^ mesons. The different corrections are described in the following. 

After this procedure, the track momentum distributions are defined as the momentum 
distributions of all promptly produced stable charged particles and those produced in the 
decays of particles with lifetimes shorter than 3x 10~^° sec, corrected for initial state radiation. 
This means that charged decay products from Kg, hyperons and weakly decaying b and c 
fiavoured hadrons are included in the definition, regardless of how far away from the interaction 
point the decay actually occurred. 



3.1 Track momentum resolution and efficiency 

The number iV°^^^"<^d of tracks in a tag sample was measured in 22 different Xp bins The 
corrected distribution for a given tag sample is 

^corrected = ^ ^ ^ {wI.N;^--'') . (2) 
j q H 

Here, Ailj is the probability that a track measured in Xp bin j originates from a true Xp bin i. 
This correction was applied to account for the migration of the tracks between different bins 
due to the track momentum resolution. The reconstruction efficiency for tracks belonging 
to a true Xp bin i is accounted for by factors ej. Differences in the slope of the Xp spectrum 
between uds, c and b quark events lead to flavour dependent migration effects and to a flavour 
dependent efficiency. Consequently, the matrix Aifj and ej are flavour dependent and have to 
be apphed to the fraction w'j f. of observed tracks created in a g = uds, c or b event. 

These weights w'j are the normalised products of the flavour dependent fragmentation 
function Fj in an Xp bin j and the fraction of events of a primary quark q in the considered 
tag sample: 

^Ik = ^ Jr^q' ■ (3) 

l^q' Ik 

The applied weights and the obtained fragmentation functions are strongly correlated. This 
was taken into account in an iterative procedure, whereby the result of the measurement was 
used to re-calculate w'j j^ and to repeat the correction procedure until the results were stable. 

^ For the measurement of the distribution, a different binning with 29 bins was used. Apart from the 
binning, there were no differences between the analysis of the Xp and the distribution, so the measurement 
of the ^p distribution is not expHcitly described in the following sections. 
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Initial values for the weighting factors were taken from Monte Carlo, but alternative initial 
values were also tried to confirm that the results did not depend on the choice of the initial 
values. 

The values for Mij and were obtained from Monte Carlo. The diagonal elements 
of the matrix, i.e., the probability that a track measured in its true Xp bin is around 80% 
in most bins, but becomes significantly smaller for high Xp bins in c and b flavoured events. 
Values for the efficiency arc typically around ^ 90% with the exception of the lowest Xp bin 
where the efficiency is about 50%. The efficiency shows only a weak flavour dependence, the 
values differ for different flavours by less than 5%. 

The corrected number of tracks A^conccted ^^^^^i tag sample was divided by the corre- 
sponding number of tagged hemispheres N^^^^ to form a fragmentation function for each 
tagged sample: 

yycorrected 

JF.fc^ . . (4) 

' ^hemi V / 



3.2 Flavour tagging and hemisphere correlations 

0*=*= mesons with high values of xd* and secondary vertices with large values for the decay 
length significance L/gl are more likely to be found in high energy jets. The hemisphere 
containing the highest energy jet also tends to have a higher charged particle multiplicity 
and a harder track momentum spectrum than the opposite hemisphere. Consequently, the 
measured fragmentation functions in samples with high values for L/ai or xd* would be 
too soft. To correct for this bias, the whole analysis was performed separately for the case 
where the tag-hemisphere contains the highest energy jet and where this is not the case. The 
unweighted average of the two results was taken at the end. 

The dependence of the track momentum spectrum on the actual value of the decay length 
significance and on the D*^ energy was also considered. Requiring a high value for L/ctl or x^* 
reduces the phase space for gluon bremsstrahlung and thus introduces a kinematic correlation 
between the hemispheres. The effect is flavour dependent, becoming more important for 
higher values of Xp and is more pronounced for the D*='=tag than for the secondary vertex 
tag. Besides this kincmatical effect, correlations also occur due to geometrical effects: in a 
typical two jet event, the jets are back to back, thus pointing into geometrically opposite 
parts of the detector. This introduces a hemisphere correlation if the detector response is not 
uniform. In addition to the kinematical and the geometrical correlations, the difference in 
the fragmentation functions in tagged events and in unselected events had to be taken into 
account. 

All these effects were accounted for by applying correction factors for each tag sample, 
flavour and Xp bin: 

^ i^^fc (generated) 

i^" (generated) ' ^ ' 

i.e., the ratio of the generated fragmentation functions in tagged events and in events where 
the tag has not been applied. These correction factors T lead to a 10% correction for high 
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Xp values opposite a hemisphere tagged by a secondary vertex and up to a 50% correction 
for high Xp values opposite a D*^ tagged hemisphere. Technically, these factors are applied as 
corrections to the purities in the fit procedure, thus taking into account the a priori unknown 
flavour composition of the tracks in a specific Xp bin. 

3.3 Background subtraction in the D*^ samples 

The measured fragmentation functions in the D*^ signal samples J-'i^^^^^ side-band 
samples J-'^f were used to determine the fragmentation functions for genuine D*^ mesons: 

p.i _ 1 / signal _ ^BG„ ^-SB \ ff^\ 

i,k ~ _ JBG-) V ^''^ ■'^ ^iM-r i^k ■) \y ) 

The background fractions f]^'^ derived from fits to the AM distribution are listed in table 
^ To take into account differences of the hemisphere correlations for events in the signal 
and those in the side-band region, correction factors Ci^k = '^ff'^'^V^ff were applied to the 
fragmentation functions of the side-band samples, i.e., they were multiplied by the ratio of the 
correction factors T for the flavour mix of the signal and the side-band sample as predicted 
by the simulation. 



4 Fits 

A simultaneous fit was performed on the fragmentation functions J-^^l^ and J^^)^ of the three 
D*^ and the three secondary vertex tagged samples to extract the flavour dependent fragmen- 
tation functions F^"^^, and F^. The fragmentation functions obtained from samples tagged 
by D*^ decays J-^f'l , corrected for detector effects and for the combinatorial background for 
each of the three xd* regions k and each of the 22 Xp bins i were described in the fit by 

^Z" = innT)Ft + {virif)Fi (7) 

where the purities P^, are given in Table ^ and the correction factors T^^ are defined in 
equation ^ Naively, the fragmentation function corrected for detector effects for each of the 
three samples tagged by secondary vertices and each of the 22 Xp bins could be described by 

= (/r^nr'^)i^"'^ + u^ii)f: + uirit)Fi (s) 

However, the fraction of c events in all three vertex tagged samples is small (see Table 
P, hence the relative uncertainty on these fractions large. If equation || would be used, the 
vertex tagged samples would dominate the fit results due to their larger statistical weight as 
compared to the D*^ tagged samples. To ensure that the D*^ samples are used to obtain 
the charm fragmentation function, the flavour fraction f^'^^ and were replaced by (1 — 
/fc)-Ruds/(-Ruds + Rc) and (1 — fjl)Rc/{Ruds + Rc) where i?uds and -Rc are the Standard Model 
values for the branching fractions Rq = r(Z° — qq)/T{7P — >■ hadrons). The decay length 
dependence of the ratio of uds to c events were accounted for by correction factors di ^ to the 
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measured fragmentation function. The fragmentation function of the secondary vertex tagged 
samples was then described by 

d^,urJf = (1 - {flTlf))Ff- + UlTlf)Ff. (9) 

where 

p puds I p pc 
^udsc _ -"-uds-Tt ~r J^cJ^ j ^^g^ 

R\ids ~l~ -^c 

The correction factors dj.fc were derived from the momentum spectrum in Monte Carlo events 
with the ratio of uds to c events taken to be the same in all vertex samples, divided by the 
unmodified momentum spectrum with a variable ratio of uds to c events. These correction 
are of the order of 1% for most of the Xp bins except for the highest Xp bins where they exceed 
10%. 

A simultaneous fit was performed to extract F^"^^, F'^ and F^. In fact, the secondary 
vertex data using equation ^ essentially fixes F^'^^^ and and then the D*^ data provide 
F^ through equation |^, allowing equation 10 to give F^^^. The fit was based on the track 



momentum spectrum of the hemisphere opposite the tag. Therefore, the results had to be 
multiplied by a factor of two to obtain the full event fragmentation functions as they are 
shown in Fig. ^ and in Table The mean values of these distributions and their statistical 
uncertainty are: 

(xp)"^' = 0.0630 ± 0.0003 
{xpf = 0.0576 ± 0.0012 
(xp)^ = 0.0529 ± 0.0001. 

The results can be compared with results for the inclusive fragmentation function which 
was obtained from the track momentum spectrum of all events without considering any flavour 
tagging. These results are shown in Table | and the mean value of the distribution was found 
to be: 

= 0.05938 ± 0.00002. 

The results for the C,p distribution are shown if Fig. ^ and Table ^. To determine the 
positions of the maxima, skewed Gaussians, i.e., combinations of two Gaussians with 
different widths to the left and to the right from the centre were fltted to these distribution 
as motivated by the next-to-leading-log (NLLA) approximation ||^. Following the procedure 
in |T^, the flt was performed in the region 2.2 < < 5.0. The results for the positions of the 
maxima with their statistical uncertainties are: 

^o"*^' = 3.74 ± 0.06 
Co' = 3.63 ± 0.16 
= 3.55 ± 0.01. 

Again, these flavour dependent results can be compared with the results of the inclusive 
Cp distribution as obtained from all events without any flavour tagging. The results are shown 
in Table ^ the position of the maximum was determined to be: 

^0*°^! = 3.656 ± 0.003 . 
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5 Systematic errors and cross checks 



The systematic uncertainties affecting the above results are due to the following sources: (1) 
uncertainties of the purities of the samples tagged by secondary vertices and (2) by D*"^ decays, 
the D*^ reconstruction (3), the hemisphere correlation (4), uncertainties inherent to the cor- 
rection procedure (5) and the track and event selection (6). These were estimated from the 
difference between the central value and the result of the repeated analysis after a cut, a purity 
or the correction procedure was modified. In each case, the largest deviation was taken as the 
systematic error. 



(1) The uncertainties on the purities in the secondary vertex event samples were estimated 
using the published results from the measurement of the charged multiplicity in b, c and 
uds events 0. There, variations of the measured multiplicity due to the uncertainties 
in the b lifetime, the fragmentation of b and c quark events, the production rates and 
the mixture of b hadrons produced as well as the decay multiplicities were studied and 
can be used to derive the uncertainties on the purities. The secondary vertex sample 
purities were then varied in the range of their uncertainty and the resulting differences 
of the results for the fragmentation function was taken as the systematic error due to 
this source. 



(2) The purities and of the D*± tag bins were taken from [jlSl. They were modified 



within their systematic errors to obtain the contribution to the systematic uncertainty 
on the fragmentation function. 

(3) To study the impact of the details of the D*^ candidate selection, the analysis was 
repeated with four sets of modified selection criteria. The M^o^'^ mass window was 
increased to 1.765 GeV/c^ < M^J^'^ < 1.965 GeV/c^; the AM window was increased to 
0.141 GeV/c^ < AM < 0.150 GeV/c^; only one of the three tracks was required to have 
a ^-chamber or jet chamber end point z measurement and the cuts based on dE/dx 
were removed. The last modification lead to a reduction of the signal to noise ratio of 
more than 25%. 

(4) The correlation between hemispheres caused by the kinematic of gluon radiation is 
corrected for by the factors T(equation |) . A good description of the energy spectrum 
and the angular distribution of jets in the Monte Carlo simulation is important for a 
reliable prediction of this effect. To estimate the effects of small discrepancies between 
data and simulation in these distribution, the analysis was repeated applying weights 
to Monte Carlo events so that the energy distribution of the most energetic jet and the 
distribution of the angle between the two most energetic jets in data and Monte Carlo 
simulation agreed. Most of the resulting event weights had values between 0.95 and 1.05. 
The difference of the results with and without weighting was taken as the systematic 
uncertainty. 

(5) To estimate systematic uncertainties due to the correction for track momentum resolu- 
tion and efficiency (Section |3.1|) and corrections applied in the fitting procedure (Section 
ffl), the following two modifications to the correction procedure were tested. First, the 
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weighting factors as defined in equation | were not re-calculated in an iterative pro- 
cedure but were based on the initial values from the Monte Carlo simulation. Secondly, 
the correction factors di^k in equation ^were omitted, i.e., the secondary vertex samples 
were not corrected for the variation of the uds to c flavour fraction. 

(6) To account for imperfections in the tracking detector simulation, results were obtained 
in six different ways with modified event and track selections and variations of track 
quantities in the Monte Carlo simulation. The cut on the angle of the thrust axis 
I cos6'thrust| > 0.8 was removed; instead of accepting all tracks, a cut on | cos6'track| < 0.7 
was applied; tracks were rejected if their 2;-coordinate at the point of closest approach 
to the event origin was larger than \zq\ > 10 cm; tracks were rejected if their momentum 
was smaller than 0.250 GeV/c; the track momenta in simulated events were modified 
by an additional smearing factor leading to a degradation of the momentum resolution 
in Monte Carlo events of 10%; the simulated track momenta were shifted by 1%. Since 
these effects are flavour independent, the uncertainty due to the track and event selection 
has been set to be the same for the inclusive and the flavour dependent fragmentation 
functions. 

The systematic uncertainties from the above groups of effects were added in quadrature 
and are shown for the flavour dependent distributions in the last column of Tables ^ and ^ 
The result for the inclusive fragmentation function was obtained without flavour tagging and 
consequently does not depend on the tagging efficiency, purity or hemisphere correlation, so 
only the last two groups of effects contributed to the systematic uncertainty shown in the last 
column of Tables § and |^. For most of the Xp range, the relative systematic error is below 
5% for uds and b events and below 10% for c events. For very high momenta {xp > 0.5), the 
systematic uncertainty becomes larger. Note that in uds and c events, the systematic and 
statistical errors are roughly equal. Detailed results are shown in Table |^ for a typical low 
momentum Xp bin (0.05 < Xp < 0.06) and a typical high momentum Xp bin (0.3 < Xp < 0.4). 

In Table ||, the systematic uncertainties on the measurement of the mean value (xp) of 
the Xp distributions are shown. Further systematic checks were done for the determination of 
the position of the maximum of the C,p distributions, ^o- Instead of evaluating the position of 
the maximum using a skewed Gaussian fit, a normal Gaussian fit as motivated by LLA WE\ 



was used. Furthermore, the fit range was modified and the skewed Gaussian was fitted to the 
measured C,p distribution in the regions 2.0 < < 5.2 and 2.4 < < 4.8. The uncertainty 
obtained from this test has been added in quadrature to the previous six contributions as 
listed in Table 

To cross-check the results, the tagging methods were modified. An alternative b tag was 
applied, based on impact parameter information rather than on decay length information: 
The third largest impact parameter of a track was taken as the tag quantity. The impact 
parameter distribution of tracks from a decay are independent from the energy of the decay- 
ing particle while the decay length of a particle is proportional to its energy. Hence, this 
simple alternative method is less affected by kinematical correlations due to gluon radiation 
but affected by other systematic effects than the standard method. The result obtained with 
the alternative b tagging method is consistent with the central values within the assigned sys- 
tematic error. Another cross-check was performed using an alternative background treatment 
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in the D*^ tagged event samples. Instead of taking the D*^ purities from W^ and subtracting 



the effect of the combinatorial background with the help of a fit to the D*^ signal, purities and 
background were taken from Monte Carlo. Again, the results obtained with this alternative 
method and the central values agreed within the systematic errors. 

Effects from the binning in the tag variable have also been cross-checked. Instead of using 
three secondary vertex bins and three D*^ bins, alternative results were obtained using two 
and four secondary vertex bins and likewise two and four D*^ bins. The deviations from the 
central value were in all cases smaller than the estimated systematic error. 

The whole analysis procedure has also been tested globally with simulated events. It was 
shown that the generated fragmentation function and the result of the unfolding procedure 
agree within the statistical uncertainty. 

By integrating the fragmentation functions, the charged multiplicity in uds, c and b events 
can be obtained. The results, 

n"ds ^ 20.25 ± 0.11 ± 0.37 

= 21.55 ± 0.37 ± 0.64 

= 23.16 ± 0.02 ± 0.45 
n''^"' = 21.16 ± 0.01 ± 0.21 , 

are in good agreement with results of direct measurements of the charged multiplicities P, ^ 



23, 24 



The average of the three flavour dependent fragmentation functions can be formed, weighted 
with the Standard Model branching fractions Rnds, Rc and R\^. This combined fragmentation 
function can be compared with the results for the inclusive fragmentation function and with 
previously published OPAL results All three results show good agreement with each 

other. 



6 Results 

The results for the flavour dependent fragmentation functions for uds, c and b events as well 
as the inclusive fragmentation function are shown in Fig. ^ and in Tables ^ and ^. The mean 
values of these distributions are: 

{xX'^'' = 0.0630 ± 0.0003 ± 0.0011 
{xp)" = 0.0576 ± 0.0012 ± 0.0016 
{xp)^ = 0.0529 ± 0.0001 ± 0.0013 
(xp)'"""' = 0.05938 ± 0.00002 ± 0.00057. 

The light quark fragmentation function is found to be harder than the b quark fragmenta- 
tion function as expected due to the cascade decays of b hadrons in b quark events with more 
particles sharing the energy. This observation is also consistent with the results of comparisons 



of gluon, uds and b jets pi], p2 . 
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In Fig. I and in Tables ^ and|^, the results are presented for the = ln(l/a;p) distribution 
which emphasises the lower momenta of the spectrum. Skewed Gaussians were fitted to these 
distributions to obtain the position of their maxima: 

^o"<i^ = 3.74 ± 0.06 ± 0.21 
^o' = 3.63 ± 0.16 ± 0.31 
^o'' = 3.55 ± 0.01 ± 0.07 
^0*°^! = 3.656 ± 0.003 ± 0.115. 

The result for the inclusive distribution is in good agreement with previous results |19|, 
whereas the position of the maxima of the flavour dependent distribution is reported here 
for the first time. Part of the systematic uncertainties cancel when the ratio of the flavour 
dependent results to the inclusive result is taken: 

^o'"^7^o'"'' = 1-023 ± 0.017 ± 0.028 
^o7^o'°=i = 0.993 ± 0.044 ± 0.082 

Another indication for a flavour dependence of the distribution is given by the differences 
of the shape of the distributions in figure ^. 

In Fig. ^, the total and the flavour dependent fragmentation functions are compared with 
results from other experiments at the same centre-of-mass energy. There is good agreement 
with [| where the D*+ K" '7T~^7r~^ decay was used as well, but a different b tagging method 
and a different correction procedure was applied. Also the quoted systematic uncertainty is 
similar in size with the exception of the high momentum region, where the uncertainty in 
is smaller than in this paper. However, a direct comparison of the systematic errors is difficult 
since the error sources dominating the systematic uncertainty in the high momentum region 
in this paper (Table |^ are not explicitly considered in ||^. Also in Fig.^, the results are shown 
to be consistent with the Jetset 7.4 expectation while the Herwig 5.9 Monte Carlo program 
||25|| p| fails to describe fully the b fragmentation function. 

The measurement of the total fragmentation function in comparison to measurements at 
lower energies [2^ and at centre-of-mass energies between 130 GeV and 161 GeV [@, is 
shown in Fig.|^. Apart from the lowest Xp region, there is good agreement with the Jetset 7.4 
prediction (solid line) despite the fact that the parameters used in Jetset 7.4 were optimised 
to describe data at the Z° resonance. In the highest Xp bin, the difference between the total 
fragmentation function measured at ^/s = 14.0 GeV and at \/s = 161 GeV is of the same 
order of magnitude as the difference between uds and b fragmentation function at ^/s = 
91.2 GeV. Observing the good agreement between data and the Jetset 7.4 prediction of the 
■\/s dependence of the total fragmentation function and of the flavour dependent fragmentation 
functions at y/s = 91.2 GeV, we can use Jetset to estimate the effect of the change of the 

''The parameter set used was the same as in |2l] for Herwig 5.8, except for the value of the cluster mass 
cutoff CLMAX which has been increased from 3.40 GeV/c^ to 3.75 GeV/c^. Alternatively, studies were done 
with the default parameter set of the Herwig Monte Carlo program. But in this case, even the inclusive 
fragmentation function fails to describe the data. 
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flavour mix to the apparent scaling violation. In Fig. ^ the Jetset 7.4 prediction is shown 
when the flavour mix at all centre-of-mass energies was forced to be the same as the flavour 
mix at ^/s = 18 GeV. Although this fixed flavour mix is very different from that at at the Z° 
peak and the flavour dependent fragmentation functions differ significantly, the changes on the 
total fragmentation function at ^/s = 91.2 GeV are of the order of only two to four percent. 
This is due to a cancellation of the effect of an increased b contribution and a decreased c 
contribution at centre-of-mass energies close to the Z° resonance. 

Comparing the positions of the maxima, of the C,p distribution, the flavour dependence 
is less pronounced than for the fragmentation function at high Xp. The values for at ^/s = 
14.0 GeV and ^/s = 161 GeV differ by almost a factor of two, while the difference between 
the flavour dependent results at y/s = 91.2 GeV is one order of magnitude smaller. Like in 
the case of the fragmentation functions, an increased b contribution and a simultaneously 
decreased c contribution leads to cancellations when comparing low energy measurements and 
results at the Z° peak. This is expected on the basis of simulations where the obtained for 
the flavour mixture at a/s = 91.2 GeV and at ^/s = 18.0 GeV differ by less than one percent. 

7 Conclusions 

Flavour dependent fragmentation functions in Z° —>■ qq events have been measured separately 
for bottom, charm and light (uds) quarks and as well as for all flavours together. These 
measurements are based on OPAL data recorded between 1990 and 1995. Event samples with 
different flavour compositions were formed using reconstructed D*^ mesons and secondary 
vertices in jets. The charged particle momentum spectrum has been studied in the event 
hemisphere opposite to the tag. A simultaneous fit was performed to extract the flavour 
dependent Xp distribution as well as the flavour dependent ^p distribution. 

The fragmentation function for b quarks is significantly softer than for uds quarks. The 
fragmentation functions are well described by the Jetset 7.4 Monte Carlo program while 
Herwig 5.9 fails to describe fully the b fragmentation function. 

For the first time, flavour dependent ^p distributions have been studied. The flavour 
dependence of the position of the maximum has been determined and was found to be small 
compared with the differences of this value at different centre-of-mass energies. 
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Table 3: Fragmentation functions of uds, c and b events. The first error is statistical, the 
second systematic. 
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± 


0.05 


± 


0.42 


0.12- 


-0.14 


28.30 


± 


0.04 


± 


0.28 


0.14- 


-0.16 


21.88 


± 


0.03 


± 


0.19 


0.16- 


-0.18 


17.10 


± 


0.03 


± 


0.16 


0.18- 


-0.20 


13.54 


± 


0.03 


± 


0.14 


0.20- 


-0.25 


9.37 


± 


0.01 


± 


0.11 


0.25- 


-0.30 


5.66 


± 


0.01 


± 


0.08 


0.30- 


-0.40 


2.89 


± 


0.01 


± 


0.05 


0.40- 


-0.50 


1.208 


± 


0.004 


± 


0.036 


0.50- 


-0.60 


0.506 


± 


0.002 


± 


0.018 


0.60- 


-0.80 


0.153 


± 


0.001 


± 


0.012 


0.80- 


-1.00 


0.0199 


± 0.0003 ± 0.0044 



Table 4: Inclusive fragmentation function. The first error is statistical, the second systematic. 
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p 


1/atot ■ daVdCp 


uds events 


c events 


b events 


0.0- 


-0.2 


0.024±0.006±0.006 


0.002±0.147±0.005 


0.0034±0.0001±0.0027 


0.2- 


-0.4 


0.114±0.003±0.011 


0.005±0.123±0.008 


0.014 ±0.001 ±0.010 


0.4- 


-0.6 


0.277±0.009±0.025 


0.07 ±0.03 ±0.06 


0.066 ± 0.002 ±0.023 


0.6- 


-0.8 


0.529±0.016±0.032 


0.20 ±0.06 ±0.07 


0.188 


± 0.004 ± 0.034 


0.8- 


-1.0 


0.86 ±0.02 ±0.05 


0.52 ± 0.08 ± 0.11 


0.40 


± 


0.01 


± 


0.05 


1.0- 


-1.2 


1.31 ±0.03 ±0.06 


0.84 ±0.12 ±0.12 


0.71 


± 


0.01 


± 


0.07 


1.2- 


-1.4 


1.76 ±0.05 ±0.07 


1.43 ± 0.16 ± 0.15 


1.15 


± 


0.01 


± 


0.07 


1.4- 


-1.6 


2.22 ± 0.06 ± 0.06 


2.10 ± 0.20 ± 0.14 


1.74 


± 


0.01 


± 


0.09 


1.6- 


-1.8 


2.70 ±0.07 ±0.06 


2.88 ±0.24 ±0.18 


2.49 


± 


0.02 


± 


0.10 


1.8- 


-2.0 


3.06 ±0.08 ±0.09 


3.89 ± 0.29 ± 0.29 


3.41 


± 


0.02 


± 


0.08 


2.0- 


-2.2 


3.76 ±0.09 ±0.11 


3.70 ±0.31 ±0.34 


4.30 


± 


0.02 


± 


0.13 


2.2- 


-2.4 


4.03 ± 0.10 ± 0.13 


4.77 ±0.35 ±0.31 


5.19 


± 


0.02 


± 


0.11 


2.4- 


-2.6 


4.48 ± 0.10 ± 0.18 


5.32 ±0.37 ±0.45 


5.98 


± 


0.02 


± 


0.15 


2.6- 


-2.8 


5.12 ± 0.11 ± 0.16 


4.83 ± 0.37 ± 0.21 


6.46 


± 


0.03 


± 


0.15 


2.8- 


-3.0 


5.22 ± 0.12 ± 0.17 


6.39 ± 0.41 ± 0.31 


6.92 


± 


0.03 


± 


0.15 


3.0- 


-3.2 


5.26 ±0.13 ±0.19 


7.89 ± 0.48 ± 0.37 


7.19 


± 


0.03 


± 


0.16 


3.2- 


3.4 


6.24 ±0.12 ±0.21 


5.48 ± 0.43 ± 0.40 


7.37 


± 


0.03 


± 


0.17 


3.4- 


-3.6 


6.02 ± 0.12 ± 0.20 


6.90 ±0.44 ±0.36 


7.49 


± 


0.03 


± 


0.17 


3.6- 


-3.8 


5.89 ±0.13 ±0.26 


7.12 ± 0.45 ± 0.74 


7.49 


± 


0.03 


± 


0.16 


3.8- 


-4.0 


6.04 ±0.12 ±0.20 


6.50 ± 0.43 ± 0.40 


7.26 


± 


0.03 


± 


0.16 


4.0- 


-4.2 


5.85 ±0.13 ±0.20 


6.24 ± 0.46 ± 0.49 


6.93 


± 


0.03 


± 


0.14 


4.2- 


-4.4 


5.58 ±0.11 ±0.14 


5.67 ±0.40 ±0.35 


6.32 


± 


0.03 


± 


0.09 


4.4- 


-4.6 


5.15 ±0.11 ±0.09 


4.84 ±0.40 ±0.16 


5.76 


± 


0.03 


± 


0.08 


4.6- 


-4.8 


4.21 ± 0.12 ± 0.24 


5.76 ±0.41 ±0.80 


4.91 


± 


0.02 


± 


0.05 


4.8- 


-5.0 


3.99 ± 0.10 ± 0.14 


3.35 ± 0.36 ± 0.43 


4.24 


± 


0.02 


± 


0.07 


5.0- 


-5.2 


2.94 ± 0.10 ± 0.15 


3.34 ± 0.33 ± 0.45 


3.30 


± 


0.02 


± 


0.06 


5.2- 


-5.4 


2.14 ±0.10 ±0.12 


3.22 ± 0.34 ± 0.40 


2.54 


± 


0.02 


± 


0.06 


5.4- 


-5.6 


1.93 ±0.08 ±0.13 


1.36 ± 0.28 ± 0.32 


1.92 


± 


0.02 


± 


0.09 


5.6- 


-5.8 


1.43 ±0.09 ±0.23 


0.78 ± 0.32 ± 0.62 


1.42 


± 


0.03 


± 


0.16 



Table 5: = ln(l/a:p) distribution of uds, c and b events. The Erst error is statistical, the 
second systematic. 
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l/fTtot • dfrVdCp 


n D- 


-0 2 


0.0153 ± 0.0003 ± 0.0035 
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-1- 
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-3 2 


6.15 


± 


0.01 


± 


0.10 


3.2- 


-3.4 


6.38 


± 


0.01 


± 


0.11 


3.4- 


3.6 


6.52 


± 


0.01 


± 


0.12 


3.6- 


-3.8 


6.48 


± 


0.01 


± 


0.10 


3.8- 


-4.0 


6.41 


± 


0.01 


± 


0.12 


4.0- 


-4.2 


6.17 


± 


0.01 


± 


0.09 


4.2- 


-4.4 


5.75 


± 


0.01 


± 


0.06 


4.4- 


-4.6 


5.24 


± 


0.01 


± 


0.04 


4.6- 


-4.8 


4.65 


± 


0.01 


± 


0.04 


4.8- 


-5.0 


3.95 


± 


0.01 


± 


0.05 


5.0- 


-5.2 


3.110 


± 


0.005 


± 


0.051 


5.2- 


-5.4 


2.423 


± 


0.004 


± 


0.055 


5.4- 


-5.6 


1.835 


± 


0.004 


± 


0.081 


5.6- 


-5.8 


1.33 


± 


0.01 


± 


0.15 



Table 6: Inclusive 
atic. 



— ln(l/xp) distribution. The first error is statistical, the second system- 
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Effect 


inclusive 


uds 


c 


b 


(1) Purities or sec. vertex samples 




( - ) 


2.2 


( 4.0) 


0.9 


( 1-5) 


1.4 


( 6.5) 


(/J r unties or U samples 




(-) 


0.2 


( 0.3) 


0.5 


(1.0) 


< 0.1 


(<o.i) 


(3) D*^ selection 




(-) 


1.2 


( 1-8) 


2.6 


( 4.2) 


< 0.1 


(0.1) 


(4) Hemisphere correlation 




(-) 


0.2 


( 0.9) 


0.4 


( 0.6) 


0.3 


(0.1) 


(5) Correction procedure 


< 0.1 


(< 0.1) 


0.8 


( 1-7) 


1.0 


( 3.3) 


1.7 


(7.6) 


(6) Track and event selection 


1.5 


(1.9) 


1.5 


( 1-9) 


1.5 


( 1-9) 


1.5 


(1.9) 


Total systematic uncertainty 


1.5 


(1.9) 


3.0 


( 5.2) 


3.4 


( 6.0) 


2.6 


(10.2) 


Total statistical uncertainty 


0.1 


(0.2) 


2.3 


( 2.6) 


6.7 


(12.5) 


0.4 


(1.0) 



Table 7: Relative systematic and statistical uncertainties in percent on the results for 0.05 < 
Xp < 0.06 (0.3 <Xp< OA). 



Effect 


inclusive 


uds 


c 


b 


(1) Purities of sec. vertex samples 




1.4 


0.8 


1.3 


(2) Purities of D*=^ samples 




0.1 


0.3 


< 0.1 


(3) D*± selection 




0.6 


2.1 


< 0.1 


(4) Hemisphere correlation 




< 0.1 


0.1 


0.2 


(5) Correction procedure 


0.04 


0.1 


1.3 


1.7 


(6) Track and event selection 


0.96 


1.0 


1.0 


1.0 


Total systematic uncertainty 


0.96 


1.8 


2.8 


2.4 


Total statistical uncertainty 


0.03 


0.6 


2.0 


0.1 



Table 8: Relative systematic and statistical uncertainties in percent on the results for the 
mean value of the Xp distribution {xp) . 



Effect 


inclusive 


uds 


c 


b 


(1) Purities of sec. vertex samples 




0.6 


0.1 


0.2 


(2) Purities of D*^ samples 




< 0.1 


0.2 


< 0.1 


(3) D*=^ selection 




0.5 


2.0 


< 0.1 


(4) Hemisphere correlation 




0.4 


0.2 


< 0.1 


(5) Correction procedure 


0.9 


0.6 


0.3 


0.5 


(6) Track and event selection 


0.6 


0.6 


0.6 


0.6 


(7) Fit range and fit type 


3.0 


5.3 


8.4 


1.9 


Total systematic uncertainty 


3.1 


5.4 


8.7 


2.1 


Total statistical uncertainty 


0.1 


1.5 


4.1 


0.3 



Table 9: Relative systematic and statistical uncertainties in percent on the results for the 
position of the maximum of the ^p distribution, ^q. 
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Figure 1: (a): The decay length significance distribution in data (symbols) and Monte Carlo 
(solid curve). The contribution from uds and from c quarks in the Monte Carlo distribution 
has been shaded. The boundaries of the three decay length significance bins used in this 
analysis: — 10 < L/aL < I, I < L/ai <5 and b < L/ai <50 are indicated by vertical 
lines, (b) to (d): The distribution of the mass difference of the D*^ candidate and D° candidate 
in the three different x d* bins. The symbols show the data while the solid lines are the results 
of the fits described in the text. 
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Figure 2: The upper plot shows the measured fragmentation functions for uds events (filled 
symbols), c events (open squares) and b events (open triangles) as well as the inclusive frag- 
mentation function (solid line). The lower plot shows the ratio of the flavour dependent 
fragmentation functions to the inclusive fragmentation function. The error bars include sta- 
tistical and systematic uncertainties. The systematic uncertainties are correlated between bins 
as well as between flavours. 
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Figure 3: — ln(l/xp) distribution for (a) all events, (b) uds events, (c) c events and (d) b 
events. The solid lines show the results of the skewed Gaussian fitted to the distributions in 
the indicated fit range and the dashed lines show the results of a normal Gaussian fit. The 
error bars include statistical and systematic uncertainties. 
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Figure 4: Comparison of the results for the inclusive fragmentation function for this analysis 
(O) with results from ALEPH (A), DELPHI (D) and MARK II (M) at = m^o ^ 
and of the flavour dependent fragmentation function with the results from DELPHI The 
error bars include statistical and systematic uncertainties. The Jetset 7.4 predictions for the 
fragmentation function are shown as full horizontal lines and the Herwig 5.9 predictions as 
dotted horizontal lines. 
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Figure 5: Comparison of the results for the inclusive fragmentation function with results at 
different lower /|2^ and higher centre-of-mass energies J^, [I^. The error bars include statistical 
and systematic uncertainties. The solid lines show the Jetset 7.4 prediction, assuming the 
centre-of-mass energy dependence of the flavour composition as predicted by the electroweak 
theory. The dotted lines show the Jetset 7.4 prediction assuming for all energies the same 
flavour mix as at ^/s = 18 GeV. The dotted line is almost entirely hidden behind the full 
line and even at -/s = 91.2 GeV, only a negligible difference between the two curves can be 
seen because the effect of an increased b contribution is compensated largely by the effect of 
a decreased c contribution. 
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